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1.  INTRODUCTION 


In  our  initial  study  of  rare  gases  (Rg),  we  examined  the  possibility  of  storing  energy  for  extended 
periods  in  the  a^r|^  excited  electronic  state  of  the  helium  dimer  (Chabalowski  ct  al.  1988).  Our 
calculations  provided  fJic  first  reliable  estimate  of  the  lifetime  of  this  state,  which  was  found  to  be  18  s 
for  the  isolated  dimer.  We  are  continuing  to  explore  the  possibilities  of  using  the  Rg  as  a  medium  for 
high-energy  density  storage, 

Apkarian  and  coworkers  (Fajardo  and  Apkarian  1986,  1988a,  1988b)  have  chosen  solids  composed 
of  Rg  atoms  as  a  starling  point  for  studying  energy  storage  via  photo-induced,  charge  transfer  (CT) 
reactions.  This  choice  was  based,  in  part,  upon  tlie  Rg  atoms  providing  "...  tlic  necessary  conceptual 
simplicity  to  afford  a  detailed  and  rigorous  description  of  tlic  funda  acnials  of  these  proces.ses."  Part  of 
tliis  "simplicity"  is  due  to  the  ab.sencc  of  complex  ground  state  interactions.  Apkaritin  et  til.  have  studied 
spectra  of  these  CT  transitions  in  liquid  xenon  (Wiedeman,  Fajardo,  and  Apkanan  1987;  Okada, 
Wiedeman,  and  Apkarian  1989;  Fajaido  et  al.  1988)  and  in  dojx^d  Rg  solids  (Fajardo  and  Apkarian  1986, 
1988a,  1988b,  1987;  Kunttu  ct  al.  1990;  Okada,  Wiedeman,  ai  d  Apkarian  1989)  (for  various  Rg  systems) 
demonstrating  that  .spatially  .separated,  long-lived,  dijx)laron  stales  are  produced  in  tlie  Xc/Cl  system  via 
the  reaction  (Fajardo  and  Apkarian  1988b); 

Xe  ^  ilCI  +  2hv  -4  Xe*(HCl)‘. 

In  addition,  llicy  have  demonstrated  Uie  Ixsing  capabilities  ofXeF  in  crystalline  argon,  creating  a  laser 
tunable  over  an  8(Fnm  range  in  tlie  visible  region  (Katz,  Feld,  and  Apkarian  1989).  Other  combination 
of  Rg  halides  are  expected  to  Ixi  likely  candidates  as  tunable  la.scrs  in  the  visible,  ultraviolet  (UV),  and 
vacuum  UV  spectral  regions  (Schwentner  and  Apkarian  1989;.  Classical  molecular  dynamics  studies  by 
Gerber  and  coworkers  (Gcrtxir,  Alimi,  and  Apkarian  1989;  Alimi,  Brokinan,  and  Gcrlwr  1989;  Alimi, 
Gcii>ci  and  Apkarian  1989)  have  shed  .some  iighl  on  [xissible  mechanisms  lor  dissociating  F2  (Alimi, 
Gcrixjr,  and  Apkarian  1990),  Clj  (Gerber,  Alimi,  and  Apkarian  1989;  Alimi,  Brokman,  and  Gerber  1989), 
and  HI  (Alimi,  Gerber,  and  Apkarian  1989)  in  Rg  .solids.  They  found  large  cjualitative  differences  in  the 
dissociation  mechanisms  as  well  as  the  dynamical  pathways  followed  by  the  dissociated  atoms  among 
the.se  tlircc  molecules. 
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Last  and  George  (LG)  (Last  and  George  1987a,  1987b,  1988;  Last,  Kim,  and  George  1987;  Last 
et  al.  1987)  have  performed  semi-cmpirica!  (i.e.,  DIM  and  DIIS)  quantum  chemical  calculations  on  various 
(Xe,  G)  and  (Xc,  HCl)  systems,  predicting  potential  energy  surfaces,  pennanent  dipole  moments, 
electronic  transition  energies,  and  electric  transition  dipole  moments.  From  tltcsc  calculations,  tliey  made 
assignments  to  the  observed  absorption  and  emission  spectra,  provided  structural  data  for  the  complexes 
in  the  gas  and  solid  states,  and  rationalized  the  observed  bandwidths.  Their  findings  are  generally  in  good 
qualitative  agreement  with  experiment.  As  LG  have  pointed  out,  there  have  been  very  few  ab  initio 
calculations  on  related  systems. 

With  respect  to  energy  storage,  the  physical  model  for  the  long-lived,  separated  charges  is  a  halide 
negative  ion  scpaiatcd  from  die  trapped,  positively  charged  hole  state  which  appears  to  be  quasi-localized 
on  Rg  molecules  of  weakly  bound  Rg  atoms.  There  is  evidence  suggesting  that  tlic  positively  charged 
Rgji(+)  chains  arc  linear  and  prefer  n=3  in  the  .solids  (Fajardo  and  Apkarian  1988b),  witii  the  Rg-Rg 
interatomic  distance  being  shorter  than  what  is  seen  in  die  neutral  Rg  matrix.  These  separated  dipolaron 
states  can  exist  so  long  as  die  negatively  charged  halide,  X(-),  remains  separated  from  die  •^g„(+)  hole. 
But  when  the  (+)  holes  migrate  to  the  X(-),  they  rapidly  form  die  cxciplcx  Xc^(+)X(-),  releasing  the 
excitoti  energy  as  fluorescence.  The  gas  phase  equilibrium  geometry  for  die  XC2CI  cxciplcx  in  the  4^1' 
state  is  predicted  to  have  a  Xc-Xc  separation  of  3.17  A  (5.99  bolir)  (Stevens  and  Krauss  1982)  and  a 
Cl-Xe  .separation  of  6.41  bohr  (Fajardo  and  Apkarian  1986). 

Tlic  cationic  .self-trapped  holes  (STHs)  can  migrate  to  the  X(-)  by  one  ol'  the  I'ollowing  two  methods; 
(!)  either  by  tunneling,  or  (2)  by  "hopping"  (Fajardo  and  Apkarian  1988b)  (due  to  a  Ihennal  perturbation) 
which  stretches  die  Xc-Xe^(+)  bond  lengdi.s  to  die  neutral  geometry,  apparently  increasing  die  energy  of 
die  Rg„(+)  to  that  of  an  excited  hole  state  and  conducting  die  hole  to  die  X(-).  At  a  constant  temperature 
of  12  K,  the  separated  dipolarons  for  Cl  in  Xe  have  been  shown  to  exist  for  longer  than  30  hr,  with  the 
ratc-dctcrmiiiing  step  to  recombination  being  die  ability  of  the  trapped  hole  to  tunnel  (Fajardo  and 
Apkarian  1988b). 

Scliwcntncr,  Fajardo,  and  Apkarian  (1989)  prc.sciiled  the  CT  process  as  the  excitation  of  a  positive 
"hole"  (analogous  to  an  clcctren  in  an  atom)  into  a  Rydberg  orbit  with  die  halide  a'nm  acting  as  a 
negatively  charged  "nucleus"  of  infinite  mass.  Ba.sed  on  diis  model,  they  calculate  the  Rydberg  term 
values  Bj,  and  the  state  radii  for  these  Rydberg  stales.  The  progre.ssion  of  predicted  B^^  energies  fit  quite 
nicely  to  the  experimentally  derived  values  obtained  from  excitation  spectra.  To  dale,  Xc  has  Ixicn  the 
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only  Rg+halidp  matrix  to  show  such  a  Rydberg  (hole)  progression  spectra  (Schwentner,  Fajardo,  and 
Apkarian  1989).  This  delocalized  Rydberg  "hole"  description  has  been  shown  to  be  applicable  for  Xe  ;  .id 
periiaps  Kr  but  not  the  lighter  Rg  atoms  (Schwentner,  Fajardo,  and  Apkarian  1989).  In  addition,  no 
evidence  has  been  leen  experimentally  to  support  a  delocalized  hole  state  for  F  in  Rg  solids.  Xe  appears 
to  play  a  special  role  among  the  Rg  in  its  ability  to  stabilize  the  charge  separation.  This  dependence  of 
tlte  charge  transfer  process  and  stabilization  on  the  halogen  and  type  of  Rg  mauix  warrants  further 
theoretical  study. 

To  obtain  an  understanding  of  the  microscopic  processes  involved  in  formation  cf  the  exciplex,  as  well 
as  the  factors  determining  the  lifetime  of  the  separated  dipolarons,  one  needs  detailed  information  about 
the  electronic  structure  and  molecular  interactions.  This  report  describes  an  ab  initio  study  of  the 
Rg/halide  interactions  of  HCl  with  Xc  atoms.  Extending  our  previous  work  on  HQ  (Adams  and 
Qiabalowski,  to  be  published),  quantum  chemical  calculations  arc  performed  on  the  ground  and  low  1>  ng 
excited  states  of  XeQ  and  HQXe,  using  effective  core  potentials  (ECPs)  to  represent  tlte  inner  shell 
electrons,  state-averaged  complete  active  space  MCSCF  (SA-^ASSCF),  and  conllguration  interaction  (Cl) 
calculations.  The  purpose  for  the  XeCl  calculations  was  to  ascertain  how  well  the  ECPs  leproduce  known 
experimcmal  molecular  properties  as  well  as  properties  calculated  by  earlier  ab  initio  studies.  In  the 
HClXe  system,  we  are  most  interested  in  determining  the  qualitative  nature  of  the  excited  electronic  states 
and  if  a  CT  state  exists  even  in  these  small  model  systems.  The  pcmianent  electric  dipole  moments  and 
transition  dipole  moments  will  also  be  predicted  as  a  function  of  HClXe  separation.  The  transition 
moments  offer  a  semi  quantitative  measure  for  the  likelihood  of  populating  each  electronic  state. 

2.  DETAILS  OF  CALCUI.ATIONS 

All  calculations  were  carried  out  in  the  Cjv  point  group.  Tlte  atomic  basis  sets  used  arc  a  combination 
of  Gaussian-type  orbitals  and  the  ECPs  of  Wadt  and  Hay  (1985).  The  ability  of  our  codes  to  liandle  ECPs 
allows  for  the  treatment  of  systems  including  such  heavy  atoms  as  Xe.  The  Guassian  basis  consists  of 
an  uncontracted  set  of  four  s-type  primitives  (Van  Duijneveldt  1971)  and  one  p-lype  polarization  function 
(ap=0.75)  on  hydrogen,  giving  [4s, Ip].  The  chlorine  basis  contains  tlircc  nonconlractcd  s-  and  p-type 
valence  atomic  obitals  (AOs)  with  exponents  optimized  for  tlieir  use  with  tlie  ECPs  (Wadt  and  Hay  1985). 
This  is  augmented  by  a  negative  ion  function  (ap=0.049)  (Dunning  and  Hay  1977)  and  a  polarization 
function  (o^j^O.SO)  (Dunning  and  Hay  1977),  as  well  as  a  set  of  noncontracted  Rydberg  functions 
{G=Qi)25,  0p=O.O2O,  0^=0.015)  (Dunning  and  Hay  1977),  for  a  total  basis  set  of  [4s,5p,2d].  Likewise, 
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the  Xc  basis  associated  with  tJte  ECPs  (Wadt  and  Hay  1985)  consists  three  s-  and  p-type  primitives 
which  remained  uncuntracted  and  were  augmented  by  a  single  polarization  function  (a(j=0.25)  (Hay  and 
Dunning  1978),  giving  [3s,3p,ld]. 

The  molec,;lar  orbitals  (MOs)  used  as  expansion  vectors  in  the  Cl  calculations  are  obtained  from  the 
SA-CASSCF  calculations  using  the  general  second-order  density  matrix  driven  MCSCF  algorithm  of 
LA;ngsficld  (1982).  These  MOs  should  be  suitable  for  describing  both  the  ground  and  excited  states.  The 
Cl  wavefunctions  arc  generated  from  the  symbolic  matrix  element,  direct-CI  method  called  ALCHEMY 
(Liu  and  Yoshimine  1981)  by  doing  all  single  ar  double  electronic  excitations  from  a  set  of  reference 
configuration  state  functions  (CSFs).  The  details  of  the  SA-CASSCF  and  Cl  approaches  will  be  given 
in  the  appropriate  sections. 

Potential  energy  curves  (PECs)  and  electronic  transition  dipole  moments  arc  then  calculated  as  a 
function  of  atomic  positions.  The  vibrational  wavefunctions  corresponding  to  tl.c  PECs  for  the  diatomic 
systems  were  determined  by  solving  the  radial  Schroedinger  equation  for  nuclear  motion  (ignoring 
rotational  motion).  Transition  probabilities  and  radiative  lifetimes  arc  predicted  based  upon  the  previously 
mentioned  treatment.  One  potential  weakness  in  these  calculatioas  is  the  neglect  of  spin-orbit  (SO) 
interactions.  But  earlier  theoretical  work  done  by  Hay  and  Dunning  (1978)  on  XeCl  showed  that  much 
useful  qualitative  and  semi-quantitative  information  can  be  obtained  without  inclusion  of  SO  effects. 

3.  Xe-G  INTERACTIONS 

3.1  Details  of  Configuration  Interaction  Calculations.  We  have  calculated  the  PECs  for  die  1,2  TT 
and  the  l,2^n  states.  In  die  C2v  point  group  tliesc  states  transform  as  X’^fA,),  n^(Bi),  and  OyfBj),  so  the 
state  symmetries  included  in  the  SA-CASSCF  arc  1,2A,,  l,2Bi,  and  1,2B2,  with  the  weighing  »  hemc 
w=(Ll,l, 1,1,1).  Due  to  the  molecular  .symmetry,  no  loss  of  generality  occurs  by  considering  only  the  x- 
component  of  the  FI  states  in  the  CIs,  with  an  appropriate  factor  of  2  applied  when  necessary.  All  the 
valence  electrons  remained  active  in  the  Cl,  and  the  entire  set  of  virtual  orbitals  was  retained.  The  MO 
fillings  in  the  two  reference  CSFs  for  the  ^Z'^(Aj)  states  arc  ( lcr^,2a^,3a^,4a\l7c^,l7ry,27c^,27Cy;  and 
(lcr,2<T^,3a\4or^,l7t^,2Jiy,2rt^,2rty),  and  for  the  ^FljjfBj)  states  arc ( lo^,2a^3a^,4o^.l7t^,lity,27tj[,27ty) 
and  ( la^,2a^,3o^,4a^,lK^,lj:y,27t^,27i:^;).  All  single  and  double  excitations  were  done  from  tlie.se  filled 
orbitals  into  die  virtual  space. 
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The  final  C'l  wavefunciions  are  then  used  to  study  the  strong  transilion‘  2*0— >1^11  and  2^1*— ♦1^1!*, 
as  well  as  die  weak  2*0  and  2^1^— tl^fl  traasitions,  reporting  lifetimes  foremissions  f:T)m  the  v'=0 

levels  This  is  believed  to  be  the  first  report  of  theoretical  predictioas  for  the  lifetimes  of  the  weak 
traasitions.  tmpincal  spin  orbit  (SO)  effects  are  then  included  to  obtain  SO-correcicd  PECs  and  electric 
dipole  traasition  moments  for  tha'e  tran.sitions. 

3.2  kesults. 

3.2.1  Neglecting  Spin-Orbit.  The  X  Cl  system  has  already  been  studied  both  theoretically  (Hay  and 
Dunning  1978;  Huber  and  Her/berg  1979)  and  experimentally  (Inouc,  Ku,  and  Setser  1984;  Griencissen, 
Xuc-Ting,  and  Komps  1981;  Vclazco  and  Setser  1975;  Brau  and  Ewing  1975).  The  current  PECs  without 
SO  effects  can  be  seen  in  Figure  1,  and  molecular  constants  are  reported  in  Table  1  along  with  values 
from  theoretical  work  by  Hay  and  Dunning  (1978)  (HD)  and  available  experimental  values.  As  can  be 
seen  in  Table  1,  our  equilibrium  bond  lengths  and  (Og’s  are  in  reasonably  good  agreement  with  HD’s 
values.  Figure  2  gives  the  electric  dipole  transition  moments,  p^,  for  the  .strong  transitions  as  a  function 
of  R(Xc-Cl).  Near  the  Re(6.4  bohr)  for  the  2  excited  state,  the  (Og  for  the  E-Z  transition  is 
approximately  4.5  times  greater  than  the  0-11  moment.  The  transition  moments  and  1  ECs  from  this  work 
are  nearly  identical  to  those  calculated  by  HD  when  SO  is  ignored. 

Table  1  also  includes  the  results  of  our  vibrational  analysis  for  the  radiative  lifetimes.  Again,  our 
predicted  lifetimes  for  the  strong  transitions,  t(2^E^->1^E'^)=6.0  ns  and  x(2^n->l^n)=74.0  ns,  arc  very 
similar  to  HDs  without  SO  effects  (i.e.,  5.6  ns  and  64.0  ns,  respectively).  In  HD’s  study,  the  SO 
corrections  had  a  significant  affect  on  these  lifetimes,  nearly  doubling  the  aforementioned  values. 

3.2.2  Including  Spin-Orbit.  It  should  be  kept  in  mind  that  with  the  inclusion  of  SO  effects,  the  term 
symbols  and  ^fl  will  only  be  approximations  to  the  correct  state  descriptions  due  to  mixing  of 
different  A  components  of  O  through  the  spin-orbit  operator.  We  have  included  SO  couplings  in  the  same 
way  as  HD,  using  the  empirically  obtained  atomic  Cl  and  Xe"^  SO  coupling  parameters.  The  inclusion 
of  SO  effects  mixes  the  2^n-2^E'^  and  l^fl-l^E^  states.  The  resulting  wavefunciions  for  the  states,  in 
HD’s  nomenclature,  arc  as  follows: 


with  0=1/2, 


Il/j(lV)  =  Ac„vH^(lVa)  I  Bc„vH'°(l2np)  =  X. 


(1) 


Tabic  1.  Molecular  Constants  and  Radiative  Lifetimes  for  XcQ  Ignoring  SO  Effects 


(0  (mv  ') 


State 

HERE 

HP* 

HERE 

HD 

3.20 

3.25 

198.0 

190.0 

2^0 

3.08 

3.14 

196.0 

188.0 

Radiative  Lifetimes  (ns) 


TRANSITION 

HERE 

HD 

EXPERirvlENT 

2V  -  l^I"^ 

6.0 

5.6 

11.1^17.0‘' 

2^n  -  i^n 

74.0 

64.0 

131.0^.53.0‘= 

2^0  -  1 V 

13.0  ps 

2^^  -  i^n 

50.0  ps 

Okada,  Wiedeman*  and  Apkarian  (19H9). 

*’  Ftjwdo  Cl  »].  (1988). 

'■  Wtdl  ind  H*y  (198.SV,  Huy  tnd  Dunning  (1978). 


With  12=-3/2, 


Iij/id^n)  -  -Bcov'»^(>^^"«)  *  Bcov‘*^<'^np), 

III,^2(22i*)  .  Ai„„T«(2Vtt)  +  B.,„T«(22np)  -  B, 

V,/2(22n)  -  -Bi„„'P“(2^S*a)  +  Bi„„‘F"(22np). 

i3/2(i^n)  -  ‘P^d^iD  =  A, 


113/2(2^0)  =  'P°(2^n)  =  c. 


0 


sr* 


where  ‘P'ci^n;  designates  the  correction  to  the  state  due  to  the  SO  interaction  with  the  l^njy2 
state.  Within  this  four-state  (in  zeroth- order)  model,  the  and  components  of  the 

n  states  have  no  other  spin  states  with  which  to  couple,  so  the  energies  of  the  and  have  been  SO 
"corrected’’  by  HD  (and  here)  according  to  the  bllowing  simple  fomiulac  used  by  HD: 

E(l3/2)  -  Ed^n)  -  (7) 

and 

E(Il3/2)  -  -  A-xe*.  (8) 

where  A^|=2y4  cm  ’  and  ,512  cm  *  arc  the  atomic  SO  coupling  constants  for  the  Cl  and  Xc^  atoms, 

respectively.  These  values  were  taken  from  HD’s  Table  II  and  assumed  to  be  independent  of  intemuclear 
separation.  The  mixing  coefficients  A  and  B,  as  well  ;ls  the  energies  for  the  0=1/2  states,  arc  obtained 
by  diagonalizing  the  2x2  matrices  as  follows: 

-v/2X, 

E(^n)+A.^  • 


The  \  is  either  (lor  the  covalent  states)  or  X-x^^  (for  the  ionic  state.s).  Tlic  X  is  defined  as  the  SO 
splitting  parameter. 


X  “ 


(9) 


Figure  3  shows  the  rEt  s  for  the  six  .states  represented  by  Equations  l-b  with  ,St)  corrections  to  the 
energies  included  as  described  previously. 

I’hc  Rp  and  (0^,  I'or  a  selection  of  SO-corrected  stales  are  given  in  Table  2.  It  can  be  seen  that  the  SO 
coneclio.,  has  had  litllc  effect  on  the  R^.  and  <0^  values  for  the  Illi^  and  IVj^  stales.  Once  again,  oui 
results  for  R^.  and  0)^  arc  seen  to  be  very  similar  to  tlio.se  of  HD  when  SO  is  included,  with  the  current 
values  lying  slightly  closer  to  experiment  than  HD’s.  For  the  llI]/2  .slate,  our  Rj.  and  tOg  values  arc  3.17  A 


t 
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Table  2.  Molecular  Constants  and  Radiative  Lifetimes  for  XeCl  Including  SO  Effects 


r,(A) 

a)g(cm’) 

Stajc  HERE 

HD* 

Exe> 

HERE 

HD 

III1/2 

3.17 

3.22 

2.94 

195.0 

'  188.0 

195 

IVi/2 

3.12 

3.18 

— 

198.0 

189.0 

— 

Transition  Moments 

and  Radiative  Lifetimes  (x) 

HERE 

HD 

Exocriment 

Transition 

i£ 

lie 

T 

X 

*''*1/2“* '‘1/2 

2.97 

8.0 

2.75 

11.0 

ILL'  27.0^ 

***3/2“**3/2 

0.98 

95.0 

0.96 

120.0 

131.0,'  53.0^ 

*^l/2“***l/2 

0.49 

168.0 

0.50 

180.0 

— 

AE(eV) 

Transition 

HERE 

HD 

Exocriment 

***l/2^*l/2 

4.39 

4.20 

4.03' 

***3/2~**3/2 

4.05 

3.76 

3.63' 

*'^l/2"***l/2 

5.32 

5.13 

— 

*  Okula,  Wiedeman,  and  Apkarian  (,19HV). 

^  Inouc,  Ku,  and  Setser  (1984). 

Fajaido  cl  al.  (1988). 

Liu  and  Yoshiininc  (1981). 

'  Wadt  and  Hay  (1985);  Hay  and  IXinning  (1978). 

f 

in  unit.?  of  I  )ubyc. 

**  T  in  units  of  ns. 

and  195  cm  ^  respectively.  The  corresponding  experimental  values  arc  2.94  A  and  195  cm  '.  HD’s 
values  are  3.22  A  and  1 88  cm  ’,  respectively.  This  .similarity  between  the  results  from  our  study  and  from 
HD’s  is  expected  due  to  the  nearly  identical  results  obtained  between  the  two  .studies  before  including  SO 
cITccts. 
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'rtic  electric  dipole  transition  moments  over  SO  slates  are  as  follows: 


<IVi/2lMel»1/2>  =  -AcovB,oN<'*^(2V)|pJ4^(l2r)>  +  BcovA,ON<'*^(22n)|p,|'P«(l2n)>.  (11) 
and 

<»3/2IMcM3/2>  =  <‘B"(22n)|pjT°(i2n)>.  (12) 

The  results  for  the  radiative  lifetimes  (including  SO)  are  also  given  in  Table  2.  The  transition  dipole 
moments  and  AEs  were  taken  in  this  study  near  the  minima  for  the  excited  states.  These  arc  R=5.80  bohr 
for  the  113^,  R=5.90  bohr  for  the  IV  and  IR=6.00  for  the  nii;2  The  lifetime,  X,  is  calculated  from 
the  definition  of  the  Einstein  cocfTicient  for  spontaneous  emission; 

A  =  .1  -  1.063  (PiO^(AE)\l0'’  .see’’.  (13) 

X  ' 

For  the  llI]/’2-^li/2  transition,  the  lifetimes  arc  x-(8  ns;l  1  ns;l  1  ns  [Inouc,  Ku,  and  ScLscr  1984],  and  27  ns 
IGrienci.s.scii,  Xuc-Ting,  and  Komps  1981])  with  the  order  being  (1)  this  study;  (2)  HD;  and 
(3)  experiments,  respectively.  Tlic  lifetimes  for  the  transition  arc  (95  n.s;12()  n.s;131  ns  (Inouc, 

Ku,  and  Sct.scr  1984|,  and  53  ns  (Gricnci.sscn,  Xuc-Ting,  and  Komps  1981]),  and  for  the 
T=(168ns;  ISOns;  no  experiment).  For  the  llli/2~**l/2  lifetime,  the  current  result  and  that  of  HD’s  arc  in 
good  agreement  with  the  more  recent  cxfjcrimcnts  of  inouc,  Ku,  Setscr  (1984).  ITic  liletimc  for  the 
Bt/2~^i3/2  transition  from  tliis  study  lies  midway  between  the  two  experimental  values  and  21%  smaller 
than  the  theoretical  value  predicted  by  HD.  Ihc  agreement  between  the  two  theoretical  studies  is  much 
better  for  the  IVjy2~>lll/2  transition  (no  experimental  data  available),  where  ourx  is  only  7%  smaller  than 
HD’s 

Tti  '  di.scrcpancy  between  the  current  ^(113/2^13, 2)  value  and  that  of  HD  is  due  in  part  to  a  significantly 
larger  AF  predicted  in  the  current  study.  Table  2  lists  the  three  transition  energies  predicted  by  the  two 
theoretical  studies,  as  well  as  the  experimental  values  for  the  llli;2~^®i/2  Iransitions.  ITiis 

.study  predicts  AE(Il3^-+l3^)=4.05  cV,  while  HD  predicts  3.76  cV,  which  is  closer  to  the  exrx'ri mental 
value  of  3.63  cV.  As  a  general  statement,  the  level  of  Cl  tlreory  used  in  the  current  study  should  be  able 
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to  generate  more  accurate  wave  functions  tlian  the  correlation  techniques  used  by  HD,  Apparently,  the 
higher-lying  state  in  this  study  did  not  receive  the  same  level  of  correlation  corrections  as  the  lower  state, 
even  though  each  wave  function  could  very  well  be  a  better  description  ol  that  respective  electronic  state 
than  what  was  produced  in  the  Cl  calculations  of  HD. 

Since  experimental  AEs  arc  usually  fairly  accurately  known,  one  commonly  used  technique  to  predict 
a  lifetime  from  theory  is  to  combine  the  theoretical  value  for  the  transition  dipole  moment  with  tlic 
experimental  AE,  thus  eliminating  the  uncertainty  in  x  due  to  the  transition  energy.  If  we  recalculate  the 
lifetimes  for  this  study  and  HD's  study  using  this  tcclmiquc,  wc  get  t=133  ns  from  this  study  and 
x=138  ns  from  HD,  using  the  electric  dipole  transition  moments  from  Table  2.  These  are  botlt  very  close 
to  the  experimental  value  of  131  ns  predicted  by  Inoue  et  al.,  with  tlic  current  study  now  being  in  better 
agreement  with  experiment  than  HD’s  value.  The  x’s  for  the  Hli/2~*li/2'  recalculated  in  this  fashion,  are 
x(this  study)=:10.4  ns  and  x(HD)=12.1  ns,  to  be  compared  with  t(exp.)=l  1.1  ns.  This  again  puts  our  value 
of  X  in  better  agreement  with  experiment  than  HD’s  value,  although  the  agreement  with  experiment  is 
quite  satisfactory  in  both  cases. 

Tlicsc  calculations  have  shown  that  using  effective  core  potentials,  along  with  modest-sized 
singlc-and  double  excitation  Cl  wave  functions,  one  can  predict  electronic  transition  probabilities  in  good 
agreement  witli  the  experimental  values,  litis  lends  credence  to  tlie  prediction  of  such  values  for 
electronic  transitions  in  which  experimentation  has  been  unable  to  determine  the  lifetimes  in  these  doped 
Rg  sy.stcms. 

4.  HCl-Xc  INTERACTIONS 


HCl  Xe  is  treated  as  a  linear  sysrem  with  the  Xc  opposite  the  hydrogen  atom.  The  atoms  arc 
coincidental  witli  the  /.-axis  of  the  global  coordinate  system.  Calculations  were  perfonned  at  Cl  Xc 


f  R(Ha  Xc)=(4 .0,5 .{),6.0,6.4,7.(),8.0. 1 0.«>, 1 2.0, 1 4 .0, 1 6.0, 1 8 .0,20.0,26.0)  bohr  while  holding 


the  H-Cl  bond  length  at  Rg=2.409  bohr  (Huber  and  Hcr/.berg  1979).  Henceforth,  R(Ha-Xe)  will  simply 
be  referred  to  as  R.  All  uniLs  arc  atomic  unle.s.s  other  /isc  stated.  The  calculations  again  u.se  effective  core 
potentials  on  Xe  and  Cl  for  all  but  the  outer  shell  electrons.  MOs  arc  obtained  at  each  point  along  the 
PEC  from  SA-CASSCF  calculations  averaging  six  electronic  states  (i.e.,  the  two  lowest  energy  roots  from 
the  A],  B|,  and  B2  IRREPs).  The  ground  .state  XCL*)  transfonns  as  the  A,  IRREP  and  tlic  two  O 
components,  (111^,11^)  as  the  (B[,B2)  IRREPs,  rcs[x;ctivcly.  A  rn)/.cii  core  which  included  the  two  lowc.sl 
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energy  MOs  of  Aj  syinmctry  was  maintained.  Ttic  numbers  of  active  MOs  per  IRREP  were  then 
(A]=4, 81=2,62=2, A2 -0),  with  tile  active  electrons  per  IRREP  being  (A]=4,Bj=4,B2=4,A2=0). 

The  use  of  EQ^s  is  analogous  to  having  a  frozen  core  of  MOs,  sti  it  must  lie  kept  in  mind  that  the 
symmetry  designations  luj,  Ihj,  etc.,  for  tlic  MOs  label  the  lowc.st  energy,  noncore  MOs  fmm  the  a^, 
bj,...IRREPs.  Witli  IK)  furtlier  freezing  of  orbitals,  these  slate-averaged  MOs  arc  u.sed  as  expansion  vectors 
in  multircfcrcncc,  singlc-and-douhle  excitation  Cl  calculations  at  each  point  along  tlic  PECs,  In  the  CIs, 
due  to  tlic  axial  symmetry  of  tlic  system,  only  the  n^(Bi)  component  of  the  fl  states  is  calculated.  Three 
slates  were  calculated  from  the  Aj  IRREP  ;md  four  slates  from  the  B,.  llirce  reference  CSFs  of  Aj 
symmetry  were  used  and  four  of  B,  symmetry.  The  reference  CSFs  are  listed  in  Table  3,  along  wiili  the 
total  number  of  CSl's  retained  [Kr  IRREP.  Since  only  singlet  states  are  treated  in  this  study,  tlic  spin 
muitiplicity  label  will  lie  dropped  from  tlie  state  labels. 

The  typical  means  for  describing  tlic  general  characteristics  of  electronic  stales  generated  in  Cl 
calculations  is  to  refer  to  tlic  nature  of  tlic  singly  occupied  MOs  in  the  main  CSI-s  Ibr  each  stale.  Tliis 
is  es[K‘cially  useful  in  determining  whether  an  clcctniiho  state  is  Rydberg,  valence,  or  some  mixture  of 
these  characters.  But  using  state-averaged  CASSCF  MOs  means  that  lliese  orbitals,  at  least  the  ones 
included  in  tlic  active  space,  are  usually  guaranteed  to  lie  .some  Rydiicrg-valciice  mixture  if  the  states  being 
averaged  are  of  both  characters,  as  it  ui)[x;ars  to  he  in  Uiis  study.  Hie  lact  that  tlie  state-averaged  MOs 
are  designed  to  simultaneously  dc.scritie  several  ttlalcs  reduces  .somewhat  their  interpretive  capabilities. 
'I'hus,  apan  from  the  (bl lowing  general  comments  regarding  the  MOs,  the  electronic  states  will  be 
discu.s.sed  primarily  in  lenns  of  their  calculated  propcitics. 

'file  Aj  states  are  described  in  itieir  main  CSI-s  by  exciting  a  single  electron  from  tlie  o(4a,)  into 
eitlier  the  .“laj  or  ha,  a*  virtual  MOs.  'Ihe  B,  states  are  dc.scrihcd  by  an  excitation  fiom  either  llie  Jt(lb[) 
bonding  or  )i(2b,)  anlitxmcling  MOs  into  cillier  the  .Saj  or  bUj  o*  MOs.  'Ihcrefoiv,  hoth  the  E(A])  and 
u.sc  the  same  <r*(*ia[,  6aj)  virtual  MOs,  which  arc  mixed  Rydticrg-valencc  in  character  for  all 
values  of  R.  At  R=2t).(),  tlie  .'ia,  becomes  essentially  a  C!  Rydlicrg  s  orbital,  and  the  Oaj  a  HCl 
a-amibonding  orbital.  In  addititm,  tlic  cT(4a))  MO,  whicli  is  doubly  occupied  in  llie  ground  state,  is  almost 
entirely  descrilicd  by  Xc  atomic  p^  orbitals  for  nio.si  R  values.  At  tlic  very  short  distances  wlicrc 
R=(4,(),.‘).()),  tliere  is  moderate  eoiitributioii  from  tlie  Cl  s  aiui  p-type  AOs,  but  the  Xe  AOs  are  still 
dominant  even  at  lliese  close  distances.  .So  txitli  the  2A,  and  3A]  elcclronic  states  ari.se  primarily  Irom 
an  electronic  cxciiation  out  ol  atomic  Xe.  All  the  R^.’s  reixirtcd  here  arc  detemiined  from  a  simple  ilirce- 
ix)inl  lit  to  a  parabola. 


Ill  nr.-^ 
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Tabic  3.  Reference  CSFs  and  Total  Number  of  CSFs  Used  in  Defining  the  Cl  Wave  Functions  for 
IlClXc 


Orbital  Occupancies  in  Reference  CSFs 

IRREPS 

la. 

2a, 

3a, 

4a, 

5a, 

6a, 

lb. 

2b, 

3b, 

Ibj 

2b, 

3b2 

A]  S'alcs 

2 

2 

2 

2 

0 

0 

2 

2 

0 

2 

2 

0 

2 

2 

2 

1 

1 

0 

2 

2 

0 

2 

2 

0 

2 

2 

2 

1 

0 

1 

2 

2 

0 

2 

2 

0 

Total  CSFs 

-  80,372 

B]  Stales 

2 

2 

2 

2 

1 

0 

2 

1 

0 

2 

2 

0 

2 

2 

2 

2 

1 

0 

1 

2 

0 

2 

2 

0 

2 

2 

2 

2 

0 

1 

2 

1 

0 

2 

2 

0 

2 

2 

2 

2 

0 

1 

1 

2 

0 

2 

2 

0 

1'olal  CSFs 

-  1 

5,344 

4.1  'I’lic  H]  Stales.  Ilic  lour  rclercnce  CSl's  (see  fable  3)  for  the  1,2,3'Ui  slates  leprcscnl  7r->a* 
excitations.  As  discussed  previously,  the  open  shell  MOs  in  the  main  CSFs  are  a  Rydbcrg-valencc  mixture 
atid  arc  more  correctly  rclerrcd  to  as  diffuse  orbitals  rather  than  pure  valence  or  Rydberg.  In  addition, 
the  bUi  MO  is  comprised  primarily  of  AOs  on  HCl  for  R=(4.0,5.(),6.()),  witJi  increasing  Xc  character  going 
from  5.0  to  8.0,  at  which  point  the  bU]  MO  is  e.ssenlially  an  even  mix  ol  UCl  ;uid  Xe  AOs.  The  Ibj  and 
2b]  MOs  are  k  Ixmding  and  anlifXKiding  in  cl.aractcr,  respectively.  For  7.0  R  12.0,  all  four  MOs 
have  significam  contributions  Irom  AOs  on  lx)lh  Cl  and  Xe,  making  the  three  11;  states  truly  .stales  of  the 
complex  and  not  isolated  to  either  UCl  or  Xc  witlriii  lliis  range.  Figure  4  contains  the  PFCs  for  the 
electronic  stales  ol  Bj  symmetry.  It  appears  that  lire  elecinrnic  states  undergo  .significant  changes  in  their 
character  anrund  R-Ki.O  Ixihr  and  again  around  20.0  Ixrlir.  Figure  5  shows  tire  ab.solutc  value  of  the 
electric  dipole  transition  moments,  p,.,  for  Uie  l,2,3Bj<— X  transitions  as  a  function  of  R(HCl-Xe).  Tliese 
loo  suggest  major  changes  taking  place  in  the  de.st;ription  of  the  electronic  slates  anrund  16,0  and 
20.0  Ixrhr. 

4.1.1  The  1 H]  vState.  The  R^.  for  this  .stale  is  predicted  to  lie  at  8.08  Ixrhr,  which  is  0.62  bohr  longer 

than  the  gnrund  state  minimum,  Figure  5  .shows  the  Iraasition  moment  for  i’Aj  changing  rapidly 

going  from  K -4.0  to  R"  6.0  Ixihr.  The  tran.sition  moment  starts  out  at  0.21  for  K=4.0,  then  apparently 

1,5 


goes  through  zero  (with  a  sign  reversal)  near  R=6.0,  and  again  rises  to  0.18  at  R=8.0.  No  attempt  was 
made  to  verify  the  changes  in  sign  for  transition  moments,  since  it  would  serve  no  purpose  in  this  study. 
Examination  of  the  Cl  coefficients  for  the  main  CSFs  in  the  IB,  and  2Bj  states  shows  them  exchanging 
character  going  from  R=5.0  to  6.0,  which  probably  accounts  for  the  rapid  change  in  the  transition  moment. 
From  R=8.0,  the  moment  slowly  rises  to  a  value  of  0.27  at  R=26.0,  with  a  slight  deviation  from  its 
monotonic  increase  occurring  around  R=18.0.  At  R=26.0,  the  transition  moment  appears  to  be 
approaching  tiic  values  reported  for  the  A(^n)<-X  transition  in  HQ,  i.c.  (HQ)=(0.37  [Adams  and 
Chabalowski,  to  be  published],  0.35-0.39  [Van  Dishoeck,  Van  Hemert,  and  Dalgamo  1982]).  Table  4 
reports  the  traasition  dipole  moments  at  selected  points  along  the  PECs,  and  Table  5  contains  the 
permanent  dipole  moments  for  an  expanded  set  of  points  on  the  PECs.  It  is  interesting  to  look  at  the 
permanent  dipole  moment  of  the  IBj  along  the  PEC.  From  Table  5  we  sec  it  holds  rather  constant  at 
values  ranging  from  approximately  2  to  3  EXibye  for  R  ^  6.0.  These  values  are  in  close  agreement  with 
the  calculated  value  of  +2.9  D  (Adams  and  Chabalowski,  to  be  published)  for  HQ  in  its  A(’n)  state. 
A  look  at  the  vertical  transition  energies,  T,  for  selected  R-values  given  in  Table  6  show  T  decreasing 
from  8.32  cV  at  R=6.0  to  8.04  cV  at  R=26.0,  again  in  reasonable  agreement  with  the  T|.=(7.90  [Adams 
and  Chabalowski,  to  be  published],  7.84  [Bettendorff,  Pcycrimhoff,  and  Buenker  1982])  cV  reported  for 
the  A('n)<— X  transition  in  HQ. 

Based  on  the  transition  energies  and  the  permanent  and  transition  dipole  moments  in  tlic  region 
R=8.0-20.0,  the  IBj  state  might  best  be  described  as  the  A(’n)  state  in  HQ  moderately  perturbed  by 
the  xenon.  For  R  <  8.0,  the  xenon  atom  begins  to  have  a  significant  affect  on  the  properties,  which  is 
reasonable  since  the  sum  of  van  dcr  Waal’s  radii  for  0(1.8  A)  and  Xc(2.2  A)  is  7.6  bohr  (4.0  A).  The 
transition  moment  shown  in  Figure  5  appears  to  be  most  strongly  affected  by  the  close  proximity  of  the 
two  species.  The  calculated  properties  for  this  state  leads  one  to  conclude  that  its  character  is  that  of  the 
valence  A('n)  excited  state  of  HQ  perturbed  by  the  Xc.  Only  at  R=26.0  docs  this  state  truly  become 
localized  to  HQ,  and  so  only  at  R=26.0  can  the  lI3j<~X  transition  be  unambiguously  assigned  to  the 
A(’n)<— X  transition  in  HQ. 

4.1.2  The  2Bj  State.  At  R=26.0,  the  2B,  is  dc.scribcd  as  a  Cl(p^)-+a*(Ha)  excitation,  with  the  AOs 
on  Xe  making  a  verj'  little  coni  bution  to  the  .singly  occupied  MOs  in  the  main  CSFs.  The  transition 
energy  given  in  Table  6  for  R=26.0  bohr  is  T=9.93  cV,  which  lies  rca.sonably  close  to  Tp=(9,65  [Adams 
and  Chabalow.ski  1993],  9.67  IBettcndorff,  Pcycrimhoff.  and  Bucr^kcr  1982])  for  the  C(*n)<-X  Rydberg 
transition  in  HQ.  The  values  reported  in  Table  4  for  the  transition  dipole  moment  in  HQ  are  ^^,-(0.57) 
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Table  4.  Absolute  Values  of  the  Electric  Dipole  Transition  Moments,  Between  the  Ground  X(’A{) 
and  Excited  States  at  Selected  R(Ha*“Xe)  Distances 


1*/ 


Statc^ 

R=6.0^ 

R=8.0 

R=12.0 

R=18.0 

R.26.0 

HCl 

2\ 

0.69 

0.68 

0.75 

0.18 

0.01 

3^Ai 

1.01 

0.34 

0.12 

0.05 

0.00 

1% 

0.04 

0.18 

0.20 

0.25 

0.27 

0.37,'^  0.35-0.39' 

0.74 

0.77 

0.55 

0.39 

0.43 

0.57, 0.48*= 

3^B, 

0.42 

0.29 

0.19 

0.08 

0.09 

^  Transition  dipole  moments  arc  in  atomic  units. 

^  Transitions  to  the  B,  states  rciircscnt  only  the  component. 

R  is  the  distance,  in  bohrs,  between  Cl  and  Xe. 

Adams  and  Chabalowski  (to  be  published),  for  R(H-C1)=2.40  bohr. 
van  Dishoeck  et  al.  for  R(H-C1)=2.41  Nihr. 


Table  5.  Permanent  Electric  Dipok:  Moments  at  Selected  R  Values  for  the  Six  Electronic  States 
in  Tliis  Study 


Pcmiaiicnt  Electric  Dipole  Moments  (Debye) 


R(bohr) 

lAi 

2Ai 

Ml 

.^1 

Ml 

Ml 

4.0 

-0.62 

+  1  55 

+0.07 

+1.09 

+6.79 

+0.10 

5.0 

-1.40 

•0.92 

+11.61 

+1.65 

+1.00 

-6.46 

6.0 

-1.53 

-0.63 

+  15.88 

+3.09 

-0.54 

-12.34 

6.4 

-1.53 

tO.04 

+17.36 

+2.99 

+0.08 

-In. 17 

7.0 

-1..52 

+  1,11 

+19.10 

+2.83 

+  1.14 

-14,17 

8.0 

1.50 

+3.15 

+21.92 

+2.15 

+3.19 

14.86 

10.0 

-1.49 

+7,25 

+27.60 

+2.35 

+7.31 

-15.72 

12.0 

-1.50 

+  11.07 

+33.20 

+2.35 

+  11.13 

-16.88 

14.0 

-1.50 

+  15.53 

+38.42 

+2.39 

+  15.57 

-17.57 

16.0 

-1.51 

+21.82 

+43.18 

+2.37 

+21.85 

-16.73 

18.0 

1.52 

+37.09 

+47.51 

+2.22 

-7.38 

+37.20 

20,0 

-1.51 

+39,89 

+51.76 

+2.29 

-10.53 

+39.96 

26.0 

-1.54 

+64.49 

+67.38 

+2.65 

-2.85 

+64,52 

NOTE:  A  positive  dipole  moment  is  defined  as  having  the  faisitivc  mid  of  the  vector  pointing  toward  the  Xe  atom. 
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Tabic  6.  Vertical  Transition  Energies,  T,  Between  the  Ground  X  (’Aj)  and  Excited  States  at  Selected 
R(HCl***Xc)  Distances 


T(eV) 


State 

R=6.(y‘ 

R=8.0‘’ 

R=18.0 

R=26.0 

HCl 

2*A] 

8.60 

8.90 

10.81 

11.47 

3'Ai 

10.81 

11.82 

13.86 

14.48 

8.32 

8.23 

8.07 

8.04 

7.90, 7.84'^ 

2'Bj 

8.75 

8.89 

9.96 

9.93 

9.65,'  9.67‘‘ 

S'Bj 

10.39 

10.70 

10.84 

14.53 

R  is  the  disuincc,  in  hohrs,  ix  iwccn  Cl  and  Xc. 

^  Cl  Xc  distance  corresponding  tcj  lowest  energy  calculated  for  die  X(1A, )  ground  state. 

*■  Adams  and  Chabalowski  (to  be  publi.shcd),  for  R(H-C1)=2.43  bohr. 

^  Bettendorff,  I’cycrimhoff,  and  Buenker  (1982),  for  R(H-CI)=2.4()  bohr. 

[Adams  and  Chabalowski,  to  be  published],  0.48  [Van  Dishocch,  Van  Hemert,  and  Dalgamo  1982]),  to 
be  compared  with  the  current  value  of  P(.(R-26.0)=0.43  atomic  tinits.  Previous  ab  initio  Cl  calculations 
(Adams  and  Chabalowski,  to  be  published)  on  HCl  predict  tlic  permanent  dipole  moment  for  the  C('n) 
state  to  be  -0.74  D,  to  be  compared  with  the  current  value  of  -2.8  at  R=26.().  Even  though  tlic  value  for 
HClXc  is  almost  three  times  that  of  HCl,  they  arc  both  relatively  small  dipole  moments  and  both  negative. 
'Ilie  2Bj  state  hxs  llie  only  negative  permanent  dipole  of  tlic  three  B,  states  at  R=26.0.  The  HCl  and  Xc 
seem  to  have  a  rather  long-range  interaction  as  can  be  seen  in  the  changing  pcmiancnt  dipole  in  Table  7, 
as  well  as  the  transition  dipole  moments  in  Figure  5,  even  as  R  incrca.scs  beyond  20  bohr  (-10.fi  A)  ,  So 
at  R-2fi.O,  the  2B|  electronic  state  docs  appear  to  be  Rydberg  in  its  extensiveness.  Tliis,  coupled  with 
the  calculated  projxtrties,  suggests  tliat,  at  R=26.0.  the  2Bj  correlates  with  the  C('n)<-X  Rydberg 
transition  in  HCl  slightly  perturbed  by  the  Xc  atom. 

It  is  ix)w  interesting  to  examine  the  rather  large  '-hanges  in  properties  as  the  Xc  approaches  the  HCl. 
In  Figure  .5,  the  transition  moment  shows  irregular  behavior  in  moving  from  R=26  to  16  bohr,  particularly 
in  the  region  R=16^  20.  From  the  PECs  in  Figure  4,  it  appears  that  the  2B]  encounters  avoided  crossings 
with  other  Bj  states  in  iliis  region,  explaining  the  peculiar  behavior  of  the  transition  moments  in  tliis 
region.  Beginning  with  p^.=0.31  atomic  units  at  R=16,  the  transition  moment  increases  as  R  decreases, 
reaching  a  maximum  of  (i^=0.78  atomic  units  at  R=7.0,  and  finally  decreasing  to  0,22  at  K=4.0,  Tltc 
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vertical  transition  energy,  T=8.89  cV,  at  R=8.0  (see  Table  6)  show  this  transition  to  be  separated  ftom  the 
lBj<— lAj  transition  by  only  0.66  cV.  This  is  considerably  less  than  the  separation  between  the  c<^-X 
and  A<-X  transitions  in  HCl,  i.e.,  AT,j=1.75  eV,  to  which  these  two  transitions  appear  to  correlate  at 
dissociation.  One  interpretation  of  this  would  be  that  the  xenon  atom  preferentially  stabilizes  the  c  state 
over  the  A  state  in  HCl.  It  is  also  worth  noting  the  deptli  of  the  lower  minimum  in  the  double-well 
potential  shown  Ibr  this  state  at  Rg=7.04  in  Figure  4.  If  we  define  the  barrier  to  dissociation  to  be 
AE  ^R=20.0’  then  AE*=1.65  cV  (37.9  kcal/mol),  which  is  a  fairly  sizable  barrier  to  dissociation.  The 

permanent  dipole  moment  is  calculated  to  be  +3.2  D  at  R=8.0,  similar  to  the  IBj  state  and  much  too  small 
to  be  considered  as  an  electron  transfer  state.  In  summary,  for  the  region  R=6.0-20.0  the  main  CSFs  for 
2Bj  show  it  to  be  a  state  composed  of  electron  density  from  both  Cl  and  Xc,  hence  having  no  counterpart 
in  the  isolated  Xc  or  HQ.  But  at  tlic  largest  R  value  treated  here,  R=26.0  bohr,  the  2Bj  state  becomes 
a  7C->a*  transition  located  primarily  on  the  HQ  and  roughly  correlates  with  tlic  c<-X  Rydberg  transition 
in  HCl  slightly  perturbed  by  the  xenon. 


4.1.3  Tlic  3B]  State.  At  R=26.0  bohr,  tliis  state  is  predicted  to  lie  14..53  eV  above  the  ground  state. 
The  earlier  tlicoretical  studies  on  HCl  (Bettcndorlf,  Peycrimhoff,  and  Buenker  1982;  Van  Dishocck, 
Van  Hcmert,  ;ind  Dalgamo  1982;  Adams  and  Qiabalowski,  to  be  published)  did  not  report  excited 
electronic  states  at  such  high  traasition  energies,  but  even  if  such  high  lying  states  had  been  reported,  the 
3B(  could  not  be  associated  with  any  HQ  stale  due  to  its  description  in  its  main  CSF.  Tliis  CSF,  witli 
Cl  coefficient  c^=().88,  represents  an  electronic  promotion  from  MO  Ibj  into  5aj,  where  Ibj  is  essentially 
the  Xe(pjj)  AO  and  Saj  is  essentially  the  Q(.s-Rydbcrg)  AO.  This  slate  then  represents  an  ionic  or 
electronic  transfer  .state  of  H  symmetry.  One  can  .subsianiialc  this  claim  by  examining  tlic  pcrmancnl 
dipole  moments  given  in  I'ablc  5,  Its  value  at  R=26.{)  is  +64.5  D,  which  obviously  makes  this  a  very 
polar  state.  Bui  if  wc  look  at  the  traasition  dipole  moment  (Table  4)  at  R=26.(),  we  see  that  there  is 
essentially  zero  probability  for  this  electron  transfer  to  occur.  Tlie  excitation  probability  docs  incrca.se  as 
R  decreases,  such  that  at  K=16.0  the  transition  moment  is  calculated  to  be  0.18  atomic  units,  and  the 


transition  energy  lia.s  dropped  to  -11  cV,  while  Ihc  permanent  dipole  moment  still  remains  very  large  but 


ol’ opposite  sign  at  -16,73  D.  This  change  in  sign  ol  tlic  permanent  dipole  moment  in  going  from  R=18.() 
to  16.0  indicates  a  radical  change  in  the  character  of  the  3B,  state,  supporting  the  earlier  sialcment  that 
the  PECs  for  the  B|  stales  show  avoided  crossings  with  one  another  in  this  region. 


For  distances  around  R-8,0,  the  state  is  still  seen  to  be  charge  transfer  in  nature  with  a  dipole  moment 
of  -14.9  D.  Il  is  intcrcsling  that  (lie  jxisilivc  end  of  the  dijxilc  vector  now  points  towards  the  hydrogen 
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atom  and  not  the  xenon.  This  would  suggest  electron  transfer  from  the  chlorine  (or  HCl)  toward  the 
xenon.  A  look  at  the  main  CSFs  for  the  3Bj  over  the  entire  PEC  sliows  that  this  state,  like  the  2B],  is 
truly  a  state  that  exists  due  to  the  complex.  The  large  and  relatively  steady  decrease  in  the  absolute 
magnitude  of  the  transition  moment  from  R=6.0  to  R=14.0  supports  the  claim  that  the  3Bj  state  does  not 
correlate  closely  with  any  state  in  the  isolated  species  over  any  part  of  the  PEC. 

4.2  The  States.  All  the  'Aj  states  calculated  here  are  predicted  to  be  bound,  although  tlic  ground 
state  lAj  is  predicted  to  have  a  barrier  to  dissociation  of  only  -0.035  eV  calculated  from  R,.=8.36  bohr. 
ITic  semi-empirical  calculations  of  Last  and  GcotTge  (1988)  predict  an  Rg(HQ — Xe)  of  7.22  bohr  with  an 
H-Q  bond  length  of  2.4 1  bohr  in  the  linear  HQXe  ground  state  complex.  The  potential  well,  however, 
is  predicted  to  very  shallow  and  broad  as  previously  noted  by  Last  and  George  (1988)  and  supported  by 
this  work  as  shown  in  Figure  6.  Based  on  the  level  of  theory  u.scd  here,  it  is  not  possible  to  conclude 
whether  or  not  the  ground  slate  is  indeed  bound.  No  efforts  were  made  to  eonect  for  the  ba.sis  set 
superposition  error  due  to  the  quiditativc  nature  of  Uiis  study.  The  calculated  D^’s  for  all  the  states  arc 
listed  in  Table  7.  Both  the  2  and  3*Aj  states  arc  described  primarily  by  electronic  promotions  from  the 
4aj  MO  of  a  symmetry  comprised  almost  entirely  of  atomic  p^  dcasity  on  the  Xc  atom  into  the  .^aj  and 
6aj  o*  MOs.  Tlicsc  arc  the  same  a*  MOs  involved  in  the  7C-->a*  excitations  describing  the  main  CSFs 
in  the  'Bj  states. 


4.2.1  The  2Aj  State.  The  dominant  CSF  in  Uic  2A,  wave  function  for  R=6.0  to  16.0  bohr  is  a 
4aj-»5a]  single  excitation.  The  square  of  the  Cl  coefficient  for  this  CSF  is  constant  at  c^=0.86  over  this 
region.  This  reprc.scnLs  an  electron  promoted  from  an  MO  comprised  essentially  of  the  Xc(p,)  atomic 
orintals  into  the  Saj  o*  MO.  Table  5  shows  that  the  permanent  electric  dipole  moment  of  the  2A,  state 
is  nearly  identical  with  the  2B(  state  for  R>6.4-lfi.O,  and  Figure  7  shows  the  two  I’ECs  to  be  essentially 
degenerate  over  tlie  region  R=8.()  to  16,0  bohr.  The  similarities  between  these  two  states  can  be 
rationalized  in  terms  of  the  sum  of  the  van  dcr  Waals  radii  for  Cl  +  Xc  which  is  predicted  to  be  7.6  bohr 
for  the  clccl-ronic  ground  state  of  HCiXc.  A.s  G  approachc.s  Xc  along  the  z,  axis,  the  electron  clouds  for 
the  two  atoms  should  begin  penetrating  one  another  around  7.6  bohr.  To  a  rough  approximation,  outside 
this  di.stancc  the  4a]  MO  can  be  tliought  oi  as  a  perturbed  Xc(p,)  atomic  orbital.  In  tlic  2Bi  state,  tire 
electronic  excitation  arises  primarily  from  Hr;  Ibj  (n^)  MO  which  is  dcscribcfl  predominandy  by  the 
Xe(p^)  atomic  orbital  for  K)8.().  Thus,  if  the  spherical  symmetry  of  tire  xenon  atom  is  only  moderately 
pcrturlrcd  by  HCl,  an  excitation  out  of  the  daj  MO  (i.c.,  Xclp,]  AO)  should  be  nearly  indistinguishable 
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R(HCI-Xe)  (au) 


Figure  7.  Potential  energy  curves  foi’  the  2Aj  anti  7.8^  electronic  states  in  linear  HCTXe. 
(See  section  4,2.1  for  a  discussion  of  these  PECs.) 
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/able  7.  Calculated  D^’s  for  HClXe® 


D^-E(4-8.0)  -  E(R=26.0),  with  R  being  the  Cl-Xc  separation  in  bohrs. 

^  Has  a  barrier  to  dissociation  of  0.02  cV  at  R=16.0  bohr,  but  predicted  to  be  unbound  by 
definition  of  given  in  footnote  a. 


from  an  excitation  out  of  the  Ibj(ti^)  MO  (i.c.,  Xelp^]  AO)  as  is  the  ease  for  the  2B,  state,  (The  reader 
is  reminded  that  both  the  A,  and  Bj  states  arc  chantclcrizcd  by  excitatioas  into  the  same  two  a*  MOs.) 
Hence,  the  FBCs  for  the  2A,  and  2Bj  begin  to  split  apart  (see  Figure  7)  for  distances  less  than  the  sum 
of  van  dcr  Waal’s  radii  which  is  near  -8  bohr. 

Between  R=5.5  and  7.0  bohr,  tlic  2A]  acts  as  if  it  is  being  preferentially  stabilized  over  the  2Bi  state 
by  the  approach  of  the  Cl  and  Xe  atoms.  The  permanent  dipole  moment  in  the  region  R=8.()~16.()  (+3.2 
and  +21,8  D  at  R=8.0  and  16.0,  respectively)  is  too  small  for  this  to  be  considered  a  charge  transfer  state, 
however,  one  might  inlcntret  it  as  a  Rydberg  state  polarized  by  the  electro  negativity  of  llic  Cl  atom.  The 
transition  moment  remains  large  and  rather  constant,  witli  values  of  pj  =0.69  atomic  units  at  R=6.0  bohr 
to  0  57  at  transition  moment  remains  both  constant  and  large  over  such  a 

distance  also  supirorts  the  interpretation  of  this  as  a  Rydberg  state.  For  values  of  R^IO.O  bohr,  the  taj(a*) 
has  a  significantly  larger  contribution  from  tlie  Xe  s-type  AOs  than  any  AOs  on  HCl.  For  thc.se  larger 
separations,  the  suite  resembles  an  Xe  p^,— >a(.s)  Rydberg  stale  ol  the  xenon  atom  strongly  perturbed  by 
the  HCl, 

In  tlie  region  R=1 8.0  20.0,  tliere  are  apparently  avoided  crossings  occurring  amongst  the  A,  excited 
.states,  having  a  signiticam  clfect  on  the  energies  (see  Figure  i ),  pcmiancnt  dipole  momenUs,  and  traasition 
dipole  moments  (sec  Figure  8).  Finally,  at  R=26.0  bohr,  ure  wave  function  is  described  by  a  CSF 
representing  an  electntn  promoted  from  the  Xe(p^)  AOs  into  tlie  Rydberg  a(s)  MO  comprised  iilmost 
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entirely  of  Cl  s-type  AOs.  The  permanent  dipole  moment  is  +64  D,  which  is  roughly  the  value  of  unit 
positive  and  negative  charges  separated  by  26  bohr.  This  now  appears  to  be  a  charge  transfer  state, 
however,  the  probability  for  reaching  this  slate  from  the  ground  state  is  practically  zero  due  to  a  vanishing 
transition  diiKtlc  moment. 

4.2.2  The  3Aj  State.  The  3'Aj  resembles  a  tme  charge  transfer  state.  This  is  substantiated  by  the 
value  of  the  permanent  dipole  moment.  At  R=8.0,  the  dipole  is  predicted  to  be  +21.9  D,  while  two  unit 
point  charges  of  opposite  sign  separated  by  8.0  bohr  would  have  a  dipole  moment  of  +20.3  D.  The 
is  predicted  to  be  5.35  bohr,  in  good  agreement  with  the  bond  length  for  the  "4i"  ionic  state  calculated 
to  be  Rg=5.52  bohr  by  Last  and  George  (1988).  It  should  be  noted  that  their  HQ  bond  length  was 
allowed  to  vary  in  their  optimization  and  was  predicted  to  be  2.99  bohr,  which  is  0.58  bohr  longer  than 
the  fixed  bond  length  of  R(HC1)=^2.41  used  here.  Tlic  vertical  tran.sition  energy  at  R=5.0  is  T=9.62  cV, 
while  Tj.=10.86  cV  calculated  with  Rg=5.35  for  the  3A|.  'Phis  can  be  compared  to  l.ast  and  George’s 
value  ofT|.=8.43  c',  taken  between  the  analogous  minimum  energy  structures.  Again  it  should  be  pointed 
out  that  our  p ict'  acd  T^  for  the  3Aj  .state  is  undoubtedly  too  high  due  to  having  fixed  the  H-Cl  bond 
length  at  2.41  bohr.  At  R=5.(),  this  .state  is  primarily  described  by  a  CSF  dcfnicd  as  a  single  excitation 
4a,-+6ai,  a"  stale.  Tliis  repre.sents  an  electronic  excitation  from  the  Xc  atomic  AOs  into  an  MO 
localized  mo.stly  on  the  HCl  and  having  ll-Q  o-antibonding  character.  This  would  suggest  a  longer  bond 
length  for  tlte  H-Cl  in  tliis  state,  con.si.stcnt  with  the  results  of  La.st  and  George  (1988). 

As  R  incrca.scs,  the  3Aj  state  soon  becomes  the  higlicst  lying  vcitical  transition  calculated  in  this 
.study,  willi  T=1 1.82  cV  at  R=8.0.  Tlte  tnuisition  energy  steadily  increases  (except  at  R=20.0,  see  Table  6) 
to  14.48  cV  at  R=26.()  bohr.  The  3Aj+-  ’x  transition  is  calculated  to  be  v(  ry  intcn.se  for  R<7.0.  The 
transition  dipole  moment  drops  rapidly  fn)m  a  value  of  11^=2.17  atomic  units  at  R=4.0  to  |Jp=().0  at  R=:10.() 
(sec  Figure  8).  Along  with  the  changing  p^,  the  character  of  the  6aj  MO  changes  from  one  arising  from 
density  on  HCl  to  one  with  contributions  from  all  three  atoms  as  R  goes  from  6  to  16  bohr.  In  the  region 
of  avoided  crossings  around  R=18-2().(),  the  transition  moment  is  seen  to  be  small  and  finally  drops  to 
essentially  zero  at  R-26.(). 

I’his  type  of  state,  at  least  in  the  region  R-6.(>-8,()  bohr,  could  serve  as  a  model  for  studying  tlic 
ch.irgc  transfer  excitation  in  solid  xenon  matrices  doped  with  HCl  as  n^poricd  by  Apkarian  and  coworkers 
(Fajardo  and  Apkarian  1986;  1987;  1988a;  1988b;  Schwentner,  Fajardo,  and  Apkarian  1989;  Wiedeman, 
Fajardo,  and  Apkarian  1987).  As  mentioned  in  the  introduction,  the  Cl-Xe  distance  has  been  estimated 
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to  be  -^.4  bohr  in  the  gas  phase  4^r  state  of  tlie  charge  transfer  cxciplcx  XC2CI.  These  calculations 
predict  a  large  probability  of  charge  transfer  in  the  3Aj  for  R=6.(>-8.0  bohr  based  upon  the  calculated 
transition  dipole  and  permanent  dipole  moments.  This  might  tlicn  be  followed  by  dissociation  of  HQ  due 
to  the  H-Cl  o*-antibonding  character  of  tlie  6aj  MO  which  places  a  major  role  in  describing  the  3Aj  state. 

5.  CONCLUSIONS 

It  has  been  shown  that  much  iasight  into  the  XeCl  interactions  can  be  achieved  through  tlie 
combination  of  effecti’  e  core  potentials  and  wave  functions  obtained  from  state  average  MCSCF  and  Cl 
calculations.  The  results  on  XeQ  compare  favorably  with  earlier  ab  initio  work  by  Hay  and  Dunning 
(1978),  but  it  is  clear  that  the  spin-orbit  corrections  now  available  to  us  will  be  indispensable  for  attaining 
quantitative  accuracy  in  predicting  transition  energies  and  transition  strengllis. 

The  results  of  the  HO-Xc  calculations  have  produced  much  insight  into  the  fundamental  nature  of  the 
excited  states  in  this  system.  Tlie  1  and  2*Bj  states  arc  predicted  to  look  like  molecular  HCl  states 
perturbed  by  tlie  Xe  atom  at  values  of  R(HQ-Xc)  less  tlian  10  bohr.  At  dissociation,  tlic.se  states  closely 
resemble  the  A(*n)  and  C('n)  states  of  HQ  plus  a  ground  state  Xc  atom.  In  contrast,  the  S'llj  stale 
shows  characteristics  associated  witli  tlie  HClXe  molecule  as  a  whole,  represented  by  an  electron 
prornotion  of  it(XcQ)-»a(HaXe). 

The  2*Aj  state  also  shows  characteristics  as.socialcd  witli  the  entire  HClXc  molecule.  But  the  3’Aj 
.state  is  by  far  the  most  interesting  and  clearly  represents  a  charge  transfer  state  witli  a  transfer  of  an 
electron  from  the  Xc  atom  to  an  antibonding  o(HQ)  MO,  ju.st  as  predicted  from  the  experiments  of 
Apkarian  and  coworkers  in  the  halogen-doixid  xenon  (and  other  Rg)  matrices. 
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City,  Stale,  Zip  Code 

7.  If  indicating  a  Oiangc  ol  Address  or  Address  CoiTcction,  plca.se  provide  the  Current  or  Correct  address 
alxwe  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  Uiis  .sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  S  l  APLE) 


